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ABSTRACT 

We investigate the constraints on basic cosmological parameters set by the first 
compact-configm-ation observations of the Very Small Array (VSA), and other cos- 
mological data sets, in the standard inflationary ACDM model. Using the weak priors 
40 < < 90 km s'^ Mpc"! and < r < 0.5, we find that the VSA and COBE-DMR 
data alone produce the constraints fitot = 1-031q;i2j ^^b^^ = 0.029iQ;gog, Vlcdmh^ = 
0.131q;q5 and Tig — 1.04j;Q;j!|g at the 68 per cent confidence level. Adding in the type 
la supernovae constraints, we additionally find fJm = 0.32lQ'Qg and JIa = 0.7lto o7- 
These constraints are consistent with those found by the BOOMERanG, DASI and 
MAXIMA experiments. We also find that, by combining all these CMB experiments 
and assuming the HST key project limits for Hq (for which the X-ray plus Sunyaev- 
Zel'dovich route gives a similar result), we obtain the tight constraints Vim — 0.28lg 
and r^A = 0.721q °3, which are consistent with, but independent of, those obtained 
using the supernovae data. 
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1 INTRODUCTION 

One of the central aims of cosmology is to determine the 
values of the fundamental cosmological parameters that de- 
scribe our Universe. A unique opportunity to achieve this 
goal is provided by the observation of anisotropies in the 
cosmic microwave background (CMB) radiation. By com- 
paring such observations with the predictions of our current 
theories of structure formation and the evolution of the Uni- 
verse, we may place constraints on the cosmological param- 
eters that appear in these models. 

The currently most favoured theoretical model for de- 



scribing our Universe is based on the idea of inflation (?), 
which provides a natural mechanism for producing initial 
density fluctuations described by a power-law spectrum with 
a slope close to unity. The simplest versions of inflation also 
predict the Universe to be spatially flat. The initial spec- 
trum of adiabatic density fluctuations is modulated through 
acoustic oscillations in the plasma phase prior to recombina- 
tion and the resulting inhomogeneities are then imprinted as 
anisotropies in the CMB. In the basic inflationary scenario, 
the CMB temperature anisotropies are predicted to follow 
a multivariate Gaussian distribution, and so may be com- 
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pletely described in terms of their angular power spectrum. 
Moreover, the acoustic oscillations in the plasma phase lead 
to a characteristic series of harmonic peaks in the power 
spectrum, which are a robust indicator of the existence of 
fluctuations on super-horizon scales. 

Although the presence of acoustic peaks in the CMB 
power spectrum is a generic prediction of inflationary mod- 
els, detailed features of the power spectrum, such as the 
relative positions and heights of the peaks, depend strongly 
on a wide range of cosmological parameters, sec e.g. ?. In- 
deed, this sensitivity to the parameters is the reason why 
observations of the CMB provide such a powerful means of 
constraining theoretical models. 

Thus measurement of the CMB power spectrum is a 
major goal of observational cosmology and numerous exper- 
iments have provided estimates of the spectrum on a range 
of angular scales. It is only recently, however, that obser- 
vations by the BOOMERanG (?), DASI (?) and MAXIMA 
(?) experiments have provided measurements of the CMB 
power spectrum with sufficient accuracy over a wide range 
of scales to allow tight constraints to be placed on a wide 
range of cosmological parameters (see, for example, ?). This 
parameter estimation process is performed by comparing the 
observed band-powers with a wide range of theoretical power 
spectra corresponding to different sets of values of the cos- 
mological parameters, which can be accurately calculated 
using the Cmbfast (?) or Camb (?) software packages. The 
comparison of the observed and predicted power spectrum 
is usually carried out in a Bayesian context by evaluating 
the likelihood function of the data as a function of the cos- 
mological parameters. 

In this paper we perform this parameter estimation pro- 
cess using, as the main CMB datasets, the flat band-power 
estimates of the CMB power spectrum measured by the Very 
Small Array (VSA) in its compact configuration, which has 
boon described in the sequence of earlier papers ?, ? and ? 
(Papers I, II and III), and the COBE-DMR band-powers for 
\aw-l normalisation. We also combine the VSA data with 
other recent CMB experiments, and constraints from the 
HST Key Project on Ho and observations of type la super- 
novae, to tighten further the constraints on the cosmological 
parameters. Two different methods are used to perform the 
parameter estimation procedure. First, we employ the tra- 
ditional technique of evaluating the likelihood function on a 
large grid in parameter space. Second, we consider a more 
flexible approach in which the likelihood function is explored 
by Markov-Chain Monte Carlo (MCMC) sampling. The lat- 
ter method has a great potential in terms of expanding the 
dimension of the parameter set which can be investigated. 
Here we use it to demonstrate the robustness of the results 
from the standard grid approach and also to provide a novel 
visualisation of the range of uncertainty in our parameter 
estimates. 



2 MODELS, PARAMETER SPACE AND 
METHODS 

In the analyses presented in this paper we restrict our atten- 
tion to the case in which the initial fluctuations are adiabatic 
with a simple power-law spectrum; such perturbations are 
naturally produced in the standard single-fleld inflationary 



model. Moreover, as is now common practice, we consider 
models in which the contents of the Universe are divided 
into three components: ordinary baryonic matter; cold dark 
matter (CDM), which interacts with baryonic matter solely 
through its gravitational effect; and an intrinsic vacuum en- 
ergy. The present-day contributions of these components, 
measured as a fraction of the critical density required to 
make the Universe spatially flat, are denoted by fib, fJcdm 
and f^A respectively. It is possible that some of the dark 
matter may, in fact, be in the form of relativistic neutrinos 
(hot dark matter) , but the presence of a hot component has 
a negligible effect on the power spectrum, given the sensitiv- 
ity and angular resolution of current CMB experiments (?). 
We therefore assume that all dark matter is cold and set 
= 0. 

Following the current theoretical expectation (?), we 
also assume that the contribution of tensor mode perturba- 
tions is very small compared with the scalar fluctuations, 
and so we ignore their effects. This assumption is consistent 
with current observations. Since tensor modes contribute 
primarily to low multipoles, i, the only existing measure- 
ment that would be particularly sensitive to their presence 
is the level of the CMB power spectrum at low-£ observed 
by the COBE-DMR experiment (?). If the tensor compo- 
nent made up a large fraction of this observed power, the 
value of the spectral index Us for scalar perturbations would 
need to exceed unity by a considerable margin in order to 
provide the level of power at higher £ measured by numer- 
ous other CMB experiments. Such a large value of Us is, 
however, ruled out by large-scale structure studies (?). Nev- 
ertheless, it must be remembered that this argument only 
holds if the initial perturbation spectrum is indeed a simple 
power- law. 

Given the assumptions outlined above, there remain 
seven degrees of freedom in the description of the stan- 
dard inflationary CDM model. The paramctcrisation of this 
seven-dimensional model space can be performed in numer- 
ous ways, but we shall adopt the most common choice, which 
is deflned by the following parameters: the physical density 
of baryonic matter {Q,f,h^ = Wb); the physical density of 
CDM (Ocdm/i^ = ij^cdm); the vacuum energy density due to 
a cosmological constant (Oa); the total density (fitot); the 
spectral index of the initial power-law spectrum of scalar 
perturbations (ris); the optical depth to the last-scattering 
surface due to reionisation (r) ; and the overall normalisation 
of the power spectrum as measured by Q = \/5C2 / (47r) and 
is quoted relative to Qcobe, as determined from the 4-year 
COBE-DMR data by ?. This choice of parameters is similar 
to that made in the analysis of the CMB band-power mea- 
surements obtained by the BOOMERanG, MAXIMA and 
DASI experiments. We note that, in this paramctcrisation, 
the reduced Hubble parameter h {= Ho km s~^ Mpc~^/100) 
is auxiliary and is given hy h = (wb + Wcdm)/(f^tot — Ha). 

In comparing the observed CMB flat band-powers mea- 
sured by the VSA with the above multidimensional model, 
we adopt a Bayesian approach based on the evaluation of 
the likelihood function for the data as a function of the 
cosmological parameters, which for brevity we denote by 
the vector 6 = (a;b, ti^cdm, f^A, fitot , J"is, t, Q/Qcobe). To set 
proper constraints on the values of these parameters, it is 
necessary to explore the seven-dimensional model space over 
a region large enough to encompass those models with sig- 
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nificant likelihood. To that end, we employ two diflerent 
techniques to explore the likelihood function: a traditional 
approach in which the likelihood function is evaluated on 
a grid of points in parameter space, and a Markov-Chain 
Monte Carlo (MCMC) technique in which a set of samples 
are drawn from the likelihood function. These techniques are 
described below. A common requirement of both methods, 
however, is the accurate evaluation of the likelihood function 
at any given point in the parameter space, and we begin by 
describing how this computation is performed. 



2.1 Evaluation of the likelihood function 

As described in Paper III, the constraints imposed by the 
VSA compact array observations on the CMB power spec- 
trum have been obtained using the Madcow analysis pack- 
age (?). This provides a one-dimensional likelihood distri- 
bution for each flat band-power Cb, conditioned on the val- 
ues of the other band-powers at the joint maximum of the 
likelihood function. In addition the Hessian matrix at the 
joint maximum is also calculated, which may be inverted to 
obtain the elements Vbb' of the covariance matrix of the 
flat band-power estimates under the assumption that the 
likelihood function is well approximated by a multivariate 
Gaussian near its peak. 

Since the VSA measures power only on small angular 
scales {i > 150), it does not constrain the amplitude and tilt 
of the power spectrum at low-^, which leads to pronounced 
degeneracies in the parameter space. Therefore, in Section 3, 
we also include in our analysis the 28 COBE-DMR band- 
powers provided in the Radpack distribution (?). Moreover, 
in Section 4, we further include CMB band-power measure- 
ments obtained by the BOOMERanG, MAXIMA and DASI 
experiments. 

In order to compare a set of observed flat band-powers 
Cb,o with our theoretical model, it is necessary to com- 
pute the corresponding predicted values Cb,p of these band- 
powers, given a particular set of parameter values 0. If Ce is 
the corresponding theoretical power spectrum for this set of 
parameter values, the predicted value of the Bth flat band- 
power is given by 

'~-B,v = 2^ -p L« 

where Ci = 1{1+1)Ci/{2t:) and Wb(^) are the window func- 
tions for the experiment under consideration. The VSA win- 
dow functions are presented in Paper III. 

In the comparison of the observed and predicted band- 
powers it is necessary to take proper account of the uncer- 
tainties in our estimates of the Cb,o- In fact, the uncertain- 
ties in band-power estimates are, in general, non-Gaussian 
and this precludes us from calculating a simple x'^-statistic. 
It is possible, however, to make a transformation from the 
flat band-powers to a related set of 'offset log-normal' vari- 
ables for which the uncertainties are Gaussian to a very good 
approximation (?). This requires the calculation of an ad- 
ditional set of quantities xb from the data, which repre- 
sent the uncertainty due to instrumental noise. These are 
straightforwardly calculated for the VSA band-powers, and 
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Figure 1. An example of true likelihood (points) and its offset 

log-normal approximation (line) for the fiat band-power in tlie 
first VSA spectral bin for the combined data from the VSAl, 
VSA2 and VSA3 fields. 



will be available on our web-site^. The corresponding quan- 
tities for the DMR and DASI band-powers are included in 
the Radpack package, and those for MAXIMA are available 
on their website. For BOOMERanG, however, the xb values 
are not yet publicly available, and so we assume simply that 
the likelihood is a multivariate Gaussian. 

The observed band-powers are then transformed as 

Cb,o = ln(CB,o + Xb), 

and similarly for the predicted band-powers. It is straight- 
forward to show that the elements of the covariance matrix 
Vbb' for the new variables are related to the covariance ma- 
trix Vbb' of the original variables by 

Vbb' = {Cb,o + XB){V)gg,{CB',o + Xb')- 

The likelihood function is then taken to be a multivariate 
Gaussian in the transformed variables, so that 

C{e) oc exp(-ix^), 
where the misfit statistic is given by 

= (Cb,o -CB,p){V)g]^,{CB',o - Cb',p)- 

b,b' 

We find that this assumed form for the likelihood function 
provides an excellent fit to the true likelihood distribution 
for the VSA band-powers. This is illustrated in Fig. 1, in 
which we plot the true likelihood and the corresponding off- 
set log-normal approximation for the flat band-power in the 
first VSA spectral bin for the combined data from the VSAl, 
VSA2 and VSA3 fields. 

We note that the use of this transformation not only 
allows us to find the best-fit point in the parameter space 
9, but also allows us to use the value of at this point as 
a measure of goodness of fit. 



http: / /www.mrao.cam.ac.uk/telescopes/vsa 
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2.2 Exploration of the parameter space: 
grid-bsised approach 

Once one has the facihty to calculate accurately the like- 
lihood function C{0) for the cosmological parameters at 
any given point in the parameter space, one must devise a 
strategy to explore the likelihood distribution throughout 
this space. 

From the above discussion, it is clear that the evalua- 
tion of the likelihood function at each point in parameter 
space requires one to calculate the theoretical power spec- 
trum corresponding to that point. The calculation of this 
Ce spectrum is performed using version 4.0 of the Cmb- 
FAST code (?), which uses a A;-space splitting technique that 
accelerates the calculation of the spectrum by using a flat 
f2tot = 1 model to compute the high-£ multipoles for mod- 
els with fitot / 1. Nevertheless, the calculation of a typical 
spectrum for a f2tot =^ 1 model still requires around 30 sec of 
CPU time on one of the processors of the Cosmos SGI Ori- 
gin 2000 computer. Even with the 16 processors available to 
us, this computational cost severely limits the total number 
of points in parameter space at which the likelihood function 
can reasonably be evaluated. We note that, in fact, two com- 
plete grids of models (and two independent pipelines for the 
cosmological parameters estimation) were set up in Tenerife 
and Cambridge. 

The traditional approach is thus to calculate the Ce 
spectra for a grid of models that is as fine as CPU-time lim- 
itations allows, while being sufficiently large to encompass 
the entire probability distribution. For our model space, we 
calculated spectra on a six-dimensional grid correspond- 
ing to the values for each parameter given in Table 1. The 
six-dimensional grid contains 419 968 spectra. Since the dif- 
ferent values of the scalar spectral index Ug can all be ob- 
tained on a single call to Cmbfast, this required 52 496 runs 
in total and took ~ 2 days of CPU time on Cosmos. The 
overall normalisation parameter Q/Qcobe need not be pre- 
computed on a grid since it produces a simple linear scaling 
of the Ce spectra. In this parameter, the likelihood function 
was calculated at 10 points in the range 0.45 to 1.44, with a 
step size of 0.11. 

The likelihood function C{0) for a given dataset 
can be evaluated at each of these grid points as de- 
scribed in section 2.1, and the location of the maximum 
determined. For each cosmological parameter, the one- 
dimensional marginalised distribution is then obtained by 
numerically integrating over the other parameters. These 
marginal distributions are then interpolated with a cubic 
spline, and determine the constraints placed on the cosmo- 
logical parameters. 

Wherever possible, we also include the calibration un- 
certainties of the CMB experiments under consideration as 
extra parameters in our analysis. The prior on such parame- 
ters is taken as a Gaussian centred on unity, with a standard 
deviation of the appropriate width. An analytic marginalisa- 
tion is performed over this parameter, using the method pro- 
posed by ?. This analytical procedure assumes, in addition 
to a Gaussian prior on the calibration parameter, that the 
likelihood function is Gaussian. Unfortunately, in neither the 
original band-powers, nor the offset log-normal variables, are 
both these functions precisely Gaussian, and so some (small) 
approximation is introduced. In this paper, the analytical 



marginalisation is performed before transforming to offset 
log-normal variables. 

2.3 Exploration of the parameter space: MCMC 
approach 

In a parameter space of large dimensionality, a natural al- 
ternative to the grid-based approach is instead to sample 
from the likelihood distribution. An efficient procedure for 
obtaining such samples is to construct a Markov chain whose 
equilibrium distribution corresponds to the likelihood func- 
tion in parameter space (see e.g. ?). Thus after propagating 
the Markov chain for a given 'burn-in' period, one obtains 
samples from the likelihood distribution, provided the chain 
has converged. 

The MCMC sampling procedure may be implemented 
most straightforwardly by using a simplified version of the 
Metropolis-Hastings algorithm. At each step n in the chain, 
the next state 0n+i is chosen by first sampling a candidate 
point 0' from some proposal distribution 7r(0). The candi- 
date point is then accepted with probability 



7r(6»„)P(6»„4 



7r(6>n+i)P(6>„) 



where V{0) is the posterior distribution for the parameters 

and is simply the product of the likelihood £{0) and and 
some prior i)(0); we take the latter to be uniform unless oth- 
erwise stated. If the candidate point is accepted, the next 
state in the chain becomes 0„+i = 0' , but if the candidate 
point is rejected, the chain does not move, and 0n+i = 0n- 
In theory, the convergence of the chain to the limiting distri- 
bution is independent of the choice of the proposal distribu- 
tion but this choice is crucial in determining both the rate 
of convergence to the limiting distribution and the efficiency 
of the subsequent sampling. An effective approach is to set 
the proposal distribution n{0) to a multivariate Gaussian, 
centred on the current point in parameter space. 

As mentioned above, the states of the chain 0„ can be 
regarded as samples from the limiting (i.e. likelihood) dis- 
tribution only after some initial burn-in period for the chain 
to reach equilibrium. The topic of convergence is still a mat- 
ter of statistical research and no definitive formula exists for 
calculating the required length of the burn-in period. Never- 
theless, several convergence diagnostics have been proposed 
(?), which may be used as a guide. 

After burn-in, the sample density is directly propor- 
tional to the likelihood distribution, and so the samples may 
be used straightforwardly to obtain estimates of the param- 
eter values and confidence limits. In particular, one may 
easily obtain one-dimensional marginalised distributions for 
each parameter separately, obviating the need for computer- 
intensive numerical integrations. Moreover, the computa- 
tional requirements of MCMC procedures are almost inde- 
pendent of dimensionality of parameter space and thus allow 
a large number of parameters to be constrained simultane- 
ously by the data. 

The strength of MCMC methods lies in the fact that 
useful parameter estimations can be achieved with consid- 
erably fewer likelihood evaluations as compared to the tradi- 
tional grid-based methods. Nevertheless, the density of sam- 
ples must be high enough so that the estimation of the un- 
derlying posterior probability distribution is not plagued by 
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Table 1. The values of tlie yjararneters in the six (hrnensional grid of models. 
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Poisson noise and is independent of the kernel density esti- 
mation method. Usually one requires on the order of several 
tens of thousands of (accepted) samples. The efficiency with 
which these samples may be obtained depends strongly on 
the shape of the posterior likelihood function and in practice 
the basic Metropolis-Hastings algorithm can be augmented 
by the use of various speed-ups. Most notably, the sampling 
efficiency is improved by considering several simultaneous 
correlated Markov chains and by specific random point gen- 
erators that attempt to follow the shape of the likelihood dis- 
tribution posterior. The main MCMC implementation used 
here is that provided by the Bayesys software (Skilling, pri- 
vate communication), which employs several enhancements 
of the basic Metropolis-Hastings algorithm and has the abil- 
ity to sample using multiple chains. 

Ideally, one would like to calculate a theoretical power 
spectrum using one of the popular packages, such as Camb 
or Cmbfast for each sample. However, this is computation- 
ally extremely expensive. As a practical alternative, one can 
use a pre-computed grid of theoretical spectra, such as that 
discussed in section 2.2, and calculate the required spec- 
trum by suitable interpolation between neighbouring grid 
points. The density of grid points in the parameter space 
must be small enough that the dominant error in the es- 
timated cosmological parameters comes from the errors in 
the measured CMB power spectrum and not from the in- 
terpolation between the grid points. We tested the accuracy 
of the grid discussed in section 2.2 by calculating the exact 
CMB power spectrum using Cmbfast for 1000 randomly 
chosen sets of parameters and comparing them with grid 
interpolation. We find that the rms error on the predicted 
band-powers resulting from the interpolation is around 4 per 
cent, which is very small as compared with the uncertainties 
in the band-powers from the current CMB experiments. 

The Bayesys sampler provides a powerful general pur- 
pose MCMC engine, which allows one to explore compli- 
cated likelihood functions of large dimensionality. We find, 
however, that the present accuracy of the CMB experiments 
results in likelihood surfaces that are relatively smooth and 
highly convex, which can be adequately explored using a 
simple MCMC sampler based solely on the straightforward 
Metropolis-Hastings algorithm. This allows the possibility 
of tailoring an MCMC algorithm to the specific problem of 
cosmological parameter estimation from CMB band-power 
measurements, and taking advantage of our prior knowledge 
concerning which parameter combinations can be calculated 
quickly using Cmbfast or Camb and which directions in 
parameter space to avoid. Such an approach has been im- 
plemented in the software package Cosmomc (?), and al- 
lows one to perform an MCMC exploration of the parame- 
ter space in which the use of a grid is bypassed altogether. 



and the exact theoretical Ce spectrum is calculated at each 
sample point. We have made use of this approach as an ad- 
ditional cross-check of our results. 

In addition to using the MCMC approach to provide 
useful checks of the parameter constraints obtained from the 
standard grid-based method, we have exploited the insensi- 
tivity of the MCMC method to the dimensionality of the 
problem by including calibration uncertainty in our numeri- 
cal analysis, rather than performing an approximate analyt- 
ical marginalisation over it, as performed in the grid-based 
approach. A new parameter a is introduced which is the ra- 
tio of the real telescope calibration to the experimental best 
estimate. A Gaussian prior is assumed for a, with its centre 
at unity and with a standard deviation corresponding to the 
calibration uncertainty of the experiment under considera- 
tion. Whenever the MCMC algorithm requires a sample for 
a given value of a, the data are rescaled accordingly. 



3 COSMOLOGICAL PARAMETER 
CONSTRAINTS FROM THE VSA 

We first consider the constraints placed on the 
values of the cosmological parameters = 
(a;b,Wcdm,nA,ntot,ns,T, Q/Qcobb) using only the VSA 
band-powers, and the low-^ normalisation provided by the 
28 COBE-DMR band-powers provided in the RADPACK 
package. The precise VSA data used were the 10 band- 
power estimates and associated covariance matrices for 
each of the three separate 3- field mosaics VSAl, VSA2 and 
VSA3. 

3.1 Grid-based approach 

Using the approach outlined in section 2.2, we calculate 
the corresponding likelihood function C{d) over the six- 
dimensional grid summarised in Table 2 and the normal- 
isation parameter Q/Qcobe- In addition, we include the 
calibration uncertainty of the VSA band-powers as an extra 
parameter in our analysis. The prior on this parameter is 
taken as a Gaussian centred on unity, with a standard devi- 
ation corresponding to the known calibration error for the 
VSA of 7 per cent in {ATf. 

After analytically marginalising over calibration un- 
certainty, the best-fit model is characterised by the 
parameter values (ctJb, ^cdm, f^A, ^tot, ns, r, Q/Qcobe) = 
(0.020,0.06,0.7,1.0,0.9,0.0,0.87), but no particular sig- 
nificance should be attached to this model. It is the 
marginalised constraints on the individual parameters that 
are most important. Nevertheless, it is of interest to deter- 
mine the goodness-of-fit for this model. At the peak, the 
value of was found to be 47.3 for the 3 x 10 VSA plus 28 
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Figure 2. Marginalised distributions for the cosmological parameters as determined from the VSA experiment, together with a normal- 
isation constraint at low-^ from COBE-DMR. The solid lines include the effect of a weak top-hat prior on h (0.4 < h < 0.9) and an 
implicit top-hat prior on r (0 < r < 0.5) resulting from the grid boundaries. The dot-dashed lines correspond to the additional prior 
r = 0. The lines are constructed using a cubic spline interpolation through the grid points, which are shown as solid circles. 



Table 2. Cosmological parameters estimated from VSA and COBE-DMR data, using several priors. It should be noted that these 
estimates include the effect of the implicit priors on the parameters resulting from the finite grid of models (see Table 1). All the 
confidence regions correspond to the 68 per cent level. 
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COBE-DMR band-powers. Assuming that a full 7 degrees of 
freedom arc lost to the fit (which is unlikely given the form 
of the theoretical power spectra), we thus have 51 remain- 
ing degrees of freedom. The value xii = 47.3 lies at the 38 
per cent point on the xii cumulative distribution function, 
which is entirely acceptable and shows that model to be a 
good fit to the data. 

For a multidimensional likelihood function calculated 
on a grid, one has already implicitly assumed top-hat priors 
on each of the parameters, corresponding to the edges of the 
grid. It should also be bourne in mind that, in principle, if 
the grid does not encompass all of the likelihood in any pa- 
rameter, then that top-hat prior becomes relevant for all of 



the parameters. In addition to the implicit prior arising from 
the grid, we may also impose explicit priors on the cosmo- 
logical parameters values, according to our existing knowl- 
edge (or prejudices) concerning their values. In particular, 
we consider combinations of the following five priors: (i) a 
weak top-hat prior on h (0.4 < h < 0.9); (ii) a strong prior 
Gaussian on h {h = 0.72 ±0.08, ?); (iii) a weak top-hat prior 
on age (10 Gyr <t <2Q Gyr); (iv) a strong prior on optical 
depth (r = 0); (v) a prior in the (0^, nA)-plane from high- 
redshift Type lA supernovae (where f2m = Ob + Ocdm) (?). 
After adopting (combinations of) these priors, we then ob- 
tain the one-dimensional marginalised distribution for each 
cosmological parameter by direct numerical integration; the 
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successive integrals over the parameter directions are evalu- 
ated in turn by first performing a cubic spline interpolation 
onto a fine regularly-space grid of points. 

An illustration of our results is shown in Fig. 2 for the 
two cases in which wo assume prior (i) above on h, with 
and without the additional prior (iv) on r. The solid lines 
represent results treating r as a free parameter, whereas the 
dot-dashed lines correspond to setting t = 0. It is clear 
from the figure that the effect of this latter prior is minimal, 
leading only to minor changes in the constraints on Qtfi^ and 
Us , which is to be expected from the well-known degeneracy 
between these two parameters (?). 

The 68 per cent confidence limits derived from these 
marginal distributions on the cosmological parameters ujb, 
Wcdm, fttot and ris are given in the first two rows of Table 2. 
The upper and lower limit in each case is defined such that 
the corresponding interval contains 68 per cent of the total 
probability, and the likelihood function evaluated at each 
limit has the same value; the quoted best-fit value is the 
mode of the corresponding marginalised distribution. Also 
listed in the table are the 68 per cent confidence limits on 
these parameters resulting from (combinations of) the dif- 
ferent priors listed above. We note that the constraints on 
these four parameters are relatively insensitive to the inclu- 
sion of increasingly stringent priors. In particular, we see 
that the constraints on Wb are all consistent with the con- 
straint Wb = 0.020 ± 0.002 (95%) (?) resulting from the 
observed primordial hydrogen-deuterium ratio and the the- 
ory of nucleosynthesis. Depending on the prior assumed, the 
preferred value of Wcdm is typically around 0.12, with an un- 
certainty of 0.05, and lends very strong support to existing 
evidence for the existence of non-baryonic dark matter. We 
also note that, for all the priors considered, the constraints 
on the total density f2tot are consistent with the Universe 
being spatially flat. Finally, and again independently of the 
particular prior assumed, the constraints on the scalar spec- 
tral index arc consistent with the scale-invariant ris ~ 1 
(Harrison-Zel'dovich) initial power spectrum, which is pre- 
ferred by standard inflationary models. 

In the flrst four rows of Table 2, the priors are unable 
to break the well-known degeneracy of CMB data in the 
(fim, fiA)-plane, and so no constraints are given. Neverthe- 
less, with the assumption of our strong prior on h and a weak 
prior on the age of the Universe, we see that one begins to 
break this degeneracy. Indeed, one flnds that the preferred 
values of flm and Qa correspond to a roughly equal parti- 
tioning of the critical density between matter (baryonic and 
dark) and vacuum energy, although the lower limit for Qm 
and the upper limit for Qa extend some way from the best- 
flt values. Finally, the assumption of the high-redshift SNIa 
constraint (?) and our weak prior on h, succeed in cleanly 
breaking the CMB degeneracy and we find the partitioning 
of the critical density between matter and vacuum energy is 
well constrained in the approximate proportions one-third 
to two-thirds. 

3.2 MCMC approach 

Using precisely the same data as analysed in the previous 
section, we also explored the seven-dimensional parameter 
space using the MCMC techniques outline in section 2.3. 
First, we used the Bayesys algorithm and, for each sam- 
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Figure 3. The marginalised distribution for tlie calibration pa- 
rameter o after analysing the data (solid line) . Also plotted is the 
original prior placed on a (dashed line). 



pie, calculated the theoretical Ci spectrum by interpolating 
from the prc-computcd grid. Since the shape of the likeli- 
hood function is very simple, it was enough to run just eight 
simultaneous Markov chains, with each walk requiring only 
a very short burn-in period. After burn-in, wc collected 3000 
samples from each chain, thus obtaining 24000 samples in 
total from the likelihood distribution. As one would hope, 
we find that all parameter constraints calculated from these 
samples are fully consistent with those obtained above for 
each set of imposed priors. We therefore do not reproduce 
them here, although they do provide a useful check on our 
earlier results. 

As mentioned in section 2.3, to illustrate the flexibility 
of the MCMC approach, we included the calibration un- 
certainty of the VSA band-power measurements as an addi- 
tional parameter in our rmmerical analysis. The prior on this 
parameter a was taken to be a Gaussian centred at unity, 
with standard deviation 0.07. In Fig. 3, we plot this prior 
distribution, together with the marginalised distribution on 
the parameter a after analysing the data. The mean of this 
distribution lies at a = 1.00 and has a standard deviation of 
0.071. Thus, we see that, within the class of models consid- 
ered, the measured CMB band-powers are consistent with 
our estimated calibration uncertainty. 

As a second illustration of the usefulness of the MCMC 
approach, we plot in Fig. 4 a novel representation of the 
sets of cosmological models consistent with the VSA plus 
COBE-DMR data, which is produced as follows. Each of the 
MCMC samples corresponds to a theoretical Ci spectrum. 
Thus from the samples it is straightforward to construct a 
one-dimensional distribution of the power at each multipole 
£. In Fig. 4, the position of the maximum of the distribu- 
tion at each £ is shown by the solid line, while the dashed 
and dot-dashed lines indicate the 68 and 95 per cent confi- 
dence limits of the distribution respectively, determined in 
the same manner as in section 3.1. This plot assumes our 
earlier weak prior on the age of the Universe. We note that 
the first peak is very well defined, and that there is also good 
evidence for the second peak. 

As mentioned in section 2.3, we also check our cosmo- 
logical parameter constraints by using the straightforward 
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Figure 4. An illustration of the MCMC sample density on the 
(/, C;)-planc. The solid line represents the maximum of the sam- 
ple density at a given £-valne, while dashed and dot-dashed lines 
correspond to 68 per cent and 95 per cent confidence limits. The 
vertical lines are VSA data points with 68 per cent confidence 
limits. Note that samples were obtained using full covariance in- 
formation, which cannot be displayed conveniently on a single 
plot. 



COSMOMC algorithm, which is optimised for the problem at 
hand and bypasses the need for a grid. Once again our re- 
sults axe fully consistent with those found above, and provide 
another useful check on our analysis. 
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4 COMBINING VSA WITH OTHER CMB 
EXPERIMENTS 

So far we have only combined the VSA data with the COBE- 
DMR experiment in order to place limits on cosmologi- 
cal parameters from the CMB. It is clear, however, that 
tighter constraints may be obtained if we additionally in- 
clude information from other CMB experiments, in particu- 
lar BOOMERanG (?), DASI (?), and MAXIMA (?). 

In Fig. 5, we begin by simply plotting the maximum- 
likelihood estimates for Obft^, Otot, ris and ficdm/i^ obtained 
by the VSA and these other recent CMB experiments, to- 
gether with the reported 68 per cent confidence intervals 
(except for MAXIMA, for which we plot as error bars 1/2 
of the reported 95 per cent confidence limits). We note that 
the confidence limits for each experiment include the effects 
of calibration and beam uncertainty, where appropriate, and 
each experiment also assumes the COBE-DMR band-powers 
and similar weak priors to those adopted in section 3.1. 

In general, the constraints on each individual param- 
eter agree within error bars. Moreover, since each of these 
experiments uses different observing techniques and has ob- 
served different regions of the sky, wc may consider each 
measurement as an independent estimate of the correspond- 
ing parameter. Assuming further that the individual likeli- 
hoods are Gaussian, they may be immediately combined to 
produce a joint constraint on each parameter, which is also 
shown in the figure. 

In the top panel, corresponding to the parameter Q,bh^ , 
we also plot the 68 per cent confidence limits arising from 
the nucleosynthesis constraint Q,bh^ = 0.020 ± 0.001 (?). 



Figure 5. Comparison of the maximum likelihood estimates and 
68 per cent confidence limits on the cosmological parameters 
0.hh?, fitot and Us from the VSA and recent other CMB experi- 
ments (see text). 



The combined measurement from all CMB experiments is 
consistent with the BBN constraint and fully supports the 
case for a low value of the primordial deuterium abundance. 
It also favours a primordial helium mass fraction of 1^ ~ 
0.246. In the second panel, we see that all the experiments 
agree with the prediction Q,tot = 1 of standard inflationary 
models. Remarkably, the combined CMB measurement has 
an error bar of only ~ 3 per cent. The combined value for ris 
also agrees at the l-o" level with the scale-invariant Harrison- 
Zel'dovich initial power spectrum, i.e. with ris = 1. Finally, 
we see that the combined value Q.^dmh'^ is tightly constrained 
to be around 5 times larger than the value for fib/i^, which is 
a strong indication for the existence of non-baryonic matter. 

In Fig. 6 we compare the VSA-|-SNIa constraints on 
f2m and , with those published by the other recent CMB 
experiments. It is important to note, however, that the pub- 
lished DASI value does not make use of any SNIa data, so 
in order to enable a proper comparison to be made, wc have 
repeated the complete grid-based parameter estimation pro- 
cedure for DASI performed by ?, using their published co- 
variance matrices and window functions, but including the 
SNIa prior from ? and a weak prior on h. Incidentally, we 
found that the parameter constraints we derived for DASI 
before including the SNIa prior were in complete agreement 
with those obtained by ?. In Fig. 6, we plot the 68 per 
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5 DISCUSSION 
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Figure 6. Comparison of tlie maximum likcliliood estimates for 
Qm and Oa from the VSA and recent CMB experiments. A prior 
from SNIa measurements is adopted (Perlmutter et al. 1999), as 
well as a weaJj prior on Hq. 



In Figs 5 and 6 we have presented results of cosmological pa- 
rameter extraction from various CMB experiments. These 
experiments use a variety of observational techniques and 
operate at a range of frequencies and have therefore widely 
different systematic effects. Nevertheless, the figures show 
remarkable agreement between different experiments. This 
may indicate the importance of the assumed weak priors 
(which are often common) and possibly the fitting of too 
many parameters given the constraining power of individ- 
ual experiments. Indeed, the reduced x is below 1 for most 
experiments (see individual papers). Nevertheless, when one 
combines experiments (VSA, DASI, BOOMERanG, MAX- 
IMA and DMR data) most cosmological parameters become 
constrained at the level ranging between 5-20 per cent, if one 
neglects the possibility of tensor modes. This accuracy rivals 
the discriminatory power of the upcoming CMB satellite ex- 
periments, such as MAP (?), with the added bonus that the 
residual systematic effects of the various experiments are 
diluted. 

We also verify that a constraint on the vacuum energy 
component of the Universe may be obtained independently 
of the Type la supernovae data, by combining the results 
from CMB experiments alone, together with the a prior on 
h from the HST key project (?). Moreover, the resulting val- 
ues for and fim are consistent with those obtained when 
supernovae data are included. It is often assumed that CMB 
experiments cannot constrain the vacuum energy, as a result 
of the well known CMB degeneracy in the (f2m, f2A)-planc, 
and the SNIa data are usually employed to break this de- 
generacy. The supernovae data are, however, still somewhat 
plagued by uncertain systematic effects (such as extinction 
along the line of sight and the evolutionary effects due to 
metallicity), although the issues are gradually being resolved 
(?). Therefore, obtaining independent consistent results on 
the value of the cosmological constant increases our confi- 
dence in both the supernovae and CMB results. 



cent confidence limits from each experiment, and we see that 
they are all in good agreement. The combined constraint is 
f^m = and Oa = 0.70i°;°t 

As noted previously by ?, by combining all the available 
CMB datasets together with a strong prior on h, it is possi- 
ble to break cleanly the CMB degeneracy in the (nm,f^A)- 
plane without using the SNIa prior. As our CMB datasets, 
we use the VSA and COBE-DMR band-powers, the RAD- 
PACK compilation for DASI, plus the published results from 
BOOMERanG and MAXIMA. By performing a simple grid- 
based parameter estimation procedure, as outlined in sec- 
tion 2.2, we find that one cam brealc the (fim, fiA)-plane 
CMB degeneracy by assuming only our weak prior on the 
age of the Universe (10 Gyr < age < 20 Gyr) and our strong 
Gaussian prior on /i (ft = 0.72 ± 0.08) from the HST key 
project. Our resulting constraints are Q,m = Q.2S^'^]^ and 
which are similar to those derived above, 
but are independent of the SNIa data. Hence, this result is 
not subject to the usual uncertainties associated with us- 
ing high-redshift type lA supernovae as standard candles. 
We stress that a very similar result would be obtained by 
instead combining the CMB data with the prior on h from 
Sunyaev-Zel'dovich and X-ray observations of clusters (?). 



6 CONCLUSIONS 

In our analysis of the newly available data from the VSA 
compact array, and other recent cosmological results, we 
have found the following. 

• Traditional grid-based methods and Markov-Chain 
Monte Carlo (MCMC) sampling techniques have been ap- 
plied to the cosmological parameter estimation problem and 
found to yield consistent results. 

• The VSA observations, when combined with the 
COBE-DMR data and a weak top-hat prior on h (0.4 < h < 
0.9) give fitot = l-03lS;l^, n^h" = 0.029iHo9, ^^cdm/i^ = 
O.lSl^;^^ and = l-04t°;JJ at the 68 per cent confidence 
level. 

• Adding in observations of type la supernovae, the CMB 
degeneracy in the (f2m, nA)-planc may be broken to yield the 
constraints 0^ = 0.321°;°^ and ^a = 0.7ll^;^^. 

• The BOOMERanG, DASI and MAXIMA experiments, 
which have different approaches and systematics, yield con- 
sistent constraints on cosmological parameters, which is 
gratifying. Combining the results of these recent CMB ex- 
periments with the VSA data, and assuming weaks priors 
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on h and the age of the Universe, gives Qtot = l-Ollg qs; 
Qbh^ = Om^taZl, ficdm/i^ = 0.12t°;°2 and = 1.00tg;gt 

• Adding in the type la supernovae constraint to the com- 
bined CMB result, the (f2m, nA)-plane degeneracy is cleanly 
broken to give Qm = 0.33iE;:;!4 and Qa = O.TOt^'ESt 

• One can equally well break the (f2m, OA)-plane degener- 
acy without the SNIa data, by assuming the strong prior on 
h from either the HST Key Project or Sunyaev-Zel'dovich 
and X-ray observations of clusters. This yields the con- 
straints = 0.281^:^7 and Qa = 0.721^3- 
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